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ABSTRACT 

We study the non-thermal emissions in a solar flare occurring on 2003 May 
29 by using RHESSI hard X-ray (HXR) and Nobeyama microwave observations. 
This flare shows several typical behaviors of the HXR and microwave emissions: 
time delay of microwave peaks relative to HXR peaks, loop-top microwave and 
footpoint HXR sources, and a harder electron energy distribution inferred from 
the microwave spectrum than from the HXR spectrum. In addition, we found 
that the time profile of the spectral index of the higher-energy ( > 100 keV) 
HXRs is similar to that of the microwaves, and is delayed from that of the 
lower-energy ( < 100 keV) HXRs. We interpret these observations in terms of an 
electron transport model called trap-plus-precipitation. We numerically solved 
the spatially-homogeneous Fokker-Planck equation to determine electron evolu- 
tion in energy and pitch-angle space. By comparing the behaviors of the HXR 
and microwave emissions predicted by the model with the observations, we dis- 
cuss the pitch-angle distribution of the electrons injected into the flare site. We 
found that the observed spectral variations can qualitatively be explained if the 
injected electrons have a pitch-angle distribution concentrated perpendicular to 
the magnetic field lines rather than isotropic distribution. 

Subject headings: acceleration of particles — Sun: flares — Sun: X-rays, gamma 
rays — Sun: radio radiation 

1. Introduction 

Observations of hard X-rays (HXRs), microwaves, and occasionally gamma-rays in solar 
flares tell us that a significant amount of non-thermal particles are produced. Among them, 



1 Department of Earth and Planetary Science, Graduate School of Science, University of Tokyo, 7-3-1, 
Hongo, Bunkyo-ku, Tokyo, 113-0033 Japan; 

2 Institute of Astronomy, University of Tokyo, 2-21-1, Osawa, Mitaka, Tokyo, 181-0015, Japan; 



-2- 



HXR and microwave observations are believed to provide the most direct information on 
electrons. Because HXRs below ~ 100 keV are emitted primarily by electrons with energy 
below several hundred keV via bremsstrahlung radiation (]Brown!ll97ll ). whereas microwaves 
above ~ 10 GHz are emitted by electrons above several hundred keV via gyrosynchrotron 



( lRamatylll969l ; lBastianlll999l ). these two sources of emission give us information on electrons 
in two different energy ranges. Therefore, a comparative study by using both HXR and 
microwave observations is useful for discussing the physics of flare non-thermal electrons 
over a wide range of energies. 



Impulsive behavior is commonly seen in both HXR and microwave lightcurves ((Kane 



19741 ). but the two emissions do not necessarily behave identically. 



Crannell et al. 



1978 



1983: 


Bai & Dennis 


1985 





Temporally, higher- 



19970 . 



Aschwanden et al. ( 19971 ) statistically analyzed the low-pass filtered HXR lightcurves for 



78 flares observed with the Compton Gamma Ray Observatory (CGRO) and find a system- 
atic increase of time delay toward higher energy. They interpreted these time delays in terms 
of electron precipitation under Coulomb collisions. 

Spatially, microwave sources do not always coincide with HXR sources. HXRs are 
typically emitted at the footpoint regions of the flare loop (|Sakaolll994r) whereas micr o waves 



are emitted mainly at the loop-top region (IMelnikov et al 



2002 



Melnikov et all (J2002J) 



suggested that only electrons with a pancake pitch-angle distribution concentrated transverse 
to the magnetic field lines can explain the observed loop-top microwave source. 



Spectrally, ISilva et al.l (120001 ) statistically studied the correlation of the HXR and mi- 
crowave spectral indices for 57 peaks of the non-thermal emission in 27 flares. They found 
that the electron energy distribution inferred from the microwave spectrum is systemati- 
cally harder than that inferred from the HXR spectrum, and suggested that the electron 
energy distribution becomes harder towards higher energy. There are three probable ex- 
planations for such spectra: (1) two (or more) different elec tron populations wit h distinct 
physical characteristics, (2) "second-s tep acceleration" (e.g.. LBai fc Ramatylll976l ). and (3) 
"trap-plus-precipitation (TPP)" (e.g.. iMelrose fc Brownlll976l ). 

Melrose fc Brown! Jl976h presented analytic solutions of th e electron energy continuit y 



equation under two conditions: strong and weak diffusion limits (IKennel &: Petschecklll966l ). 
In the strong diffusion limit, electrons injected into the loop undergo significant scattering 
and then are quickly isotropized during the loop transit. They can escape from the loop with 
a precipitation rate proportional to their velocity, u p oc v. In the weak diffusion limit, on the 
other hand, electrons are less scattered during the transit. When the loss cone distribution 
is formed, the pitch-angle diffusion time Td, which is longer than the transit time, controls 
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the electron precipitation, yielding u p oc l/ra- The precipitation rate and the evolution of 
electrons vary, depending on which condition applies. 

There have been m any observations t hat ca n be explained in terms of the TPP model 



[e.g., ILee et al.l 120001). lYokovama et al.l (120021 ) reported the Nobeyama Radioheliograph 
(NoRH; lNakajima et al.lll994l ) observation of a flare occurring on 1999 August 28. The NoRH 
observation showed clear flare-loop structure and propagating features along the loop. They 
showed that the microwave spectrum in the optically-thin regime is hard (with spectral index 
~ 1.5) around the loop-top, and then becomes softer (spectral index ~ 3.5) toward the foot- 
points. Their observation indicates that the higher-energy electrons are efficiently trapped 
within the loop, supporting the TPP model. For this event, however, no HXR observation 
was available for comparison with the microwave observation. ISilva et al.l (120001 ) pointed out 
that the discrepancy of the energy distribution between the HXR and microwave emitting 
electrons found in their study could be explained by the TPP model. In their study, there 
was no imaging observation to confirm their suggestion. If the HXR and microwave sources 
do not coincide spatially, the discrepancy of the energy distribution between the HXR and 
microwave emitting electrons can be explained by the different spatial distribution between 
the HXR and microwave emitting electrons as a result of the TPP model. Imaging as well 
as spectral data at both HXR and microwave wavelengths are essential to confirm the role 
of TPP on the parent electrons 

In this paper we analyze the non-thermal emissions of the 2003 May 29 flare by using 
the Reuven Ramaty High Energy Solar Spectroscopic Imager (R HESSI; iLin et al.ll2002l ). the 
Nobeyama Radio Polarimeters (NoRP, iNakajima et al.l Il985l . and references theirin) and 
NoRH. RHESSI has superior spectroscopic ability from ~ 3 keV to ~ 17 MeV, providing 
the HXR spectrum from ~ 3 keV to ~ 300 keV with a spectral resolution of ~ 1 keV and 
arbitrary energy bands. In previous studies, the temporal evolution of the (HXR) spectrum 
has been considered in less detail, probably due to instrumental limitations. However, the 
temporally-resolved analysis of the spectrum is important because non-thermal emissions 
and thus non-thermal electrons are the most "time-varying" objects in solar flares. RHESSI 
enables us to analyze an accurate, temporally-resolved HXR spectrum above ~ 100 keV. 
Because the HXRs above ~ 10 keV are mainly emitted by electrons above ~ 200 keV 
( lAschwanden fc Schwartz! Il996l ). RHESSI's well-resolved spectral data below ~ 300 keV 
provides us more accurate information on electrons from tens to hundreds of keV than 
before. Combining the RHESSI HXR and NoRH/NoRP microwave spectral data allows us 
to fully cover the electrons from tens to thousands of keV. 



For a physical interpretation of the observations, we use a numerical model of TPP 
which treats the pitch-angle diffusion more generally than the analytic solutions developed 



-4- 



for the weak and strong diffusion limits. iLee fc Garyl (120001 ) performed a similar treatment 
of the electron transport to explain their microwave observation of a flare on 1993 June 3. 
We also predict the microwave and HXR emissions from the calculated electron distribution. 
Comparing these model results with the observations, we discuss electron injection and 
transport, and address how the pitch-angle distribution of the injected electrons affects the 
evolution of the trapped and precipitating electrons, and their resultant emissions. 

The paper proceeds as follows. In § [2] we present a comparative study of the non- 
thermal emissions of a solar flare occurring on 2003 May 29, by using the RHESSI HXR 
and Nobeyama microwave observations. Temporally-resolved spectra of the HXRs and mi- 
crowaves are analyzed in detail. We discuss energy-dependent delays of the time profiles of 
the spectral indices, which have not been discussed in previous studies. In §[3] we present our 
treatment of the TPP model. We numerically solve the spatially-ho mogeneous Fokker-Planck 
equation (jMcClementa Il990bl ) with the Coulomb interaction (e.g.. lLeach fc Petrosianl 1 19811 1 
and a time-dependent injection. In § H] we describe the time evolution of the trapped and 
precipitating electron distribution and the predicted microwave and HXR emissions. The 
behavior of the HXR and microwave emissions predicted by the model are compared with 
the observations, allowing us to give some constraints on the properties of flare non-thermal 
electrons. In § [5] we conclude our study. 



2. Observations 

We studied a solar flare that occurred on 2003 May 29. The GOES soft X-ray (SXR) 
level was XI. 2 (Fig. [TJ upper). This flare occurred at S07W31 at 00:50 UT, and lasted 
about 1 hour. RHESSI detected a significant amount of high energy (< 300 keV) HXRs 
during 01:00 - 01:05 UT; we define this period as the impulsive phase. NoRP and NoRH 
also observed this flare. 



2.1. Lightcurves 

Figure [T] shows the lightcurves of SXRs (upper), HXRs and microwaves (lower) during 
the impulsive phase (01:00 - 01:05 UT). The upper part of the lower plot shows the HXR 
lightcurves taken with RHESSI in three energy bands: 50-70 keV, 70-100 keV, and 100-200 
keV. The dashed lines denote the peak times of each spike in the 50-70 keV band. The HXR 
lightcurves consist of four spikes, each with durations of ~ 1 minute. The peaks of the HXR 
lightcurves in the three energy bands coincide with each other within a temporal resolution 



- 5 - 



of 2 s. 

The microwave lightcurves at 17 GHz and 35 GHz presented in the lower part of the lower 
plot in Figure [T] also show four peaks, simultaneous with those in the HXRs. This indicates 
that the 17 GHz and 35 GHz microwaves are non-thermal gyrosynchrotron emissions. We 
confirm this assertion from spectral analysis in § 12. 3[ Note that the microwave peaks are 
delayed from the HXR peaks by about 4 s. 



2.2. Images 

Figure [2] shows the spatial distribution of the emissions around 01:04:30 UT. The 
TRACE 195 A image (top left), the SOHO/MDI magnetogram (top right), and the NoRH 
34 GHz brightness temperature image (bottom left) are overlaid by the RHESSI 10-20 keV, 
50-100 keV, and 100-200 ke V contours, respecti vely. The RHESSI images are reconstructed 



with the PIXON algorithm (IHurford et al.ll2002l ) using collimators 3-9 with an accumulation 



time of 60 s (01:04:00 - 01:05:00 UT). The map of the degree of polarization at 17 GHz 
(bottom right) is overlaid on the NoRH 34 GHz brightness temperature image. 

The TRACE 195 A image (Fig. [21 top left) shows a typical two-ribbon and arcade 
structure. The HXR (10-20 keV) source is co-spatial with the brightest region in 195 A, 
indicating that the 10-20 keV HXR emission is thermal bremsstrahlung from the coronal 
plasma in the flare loop. 

The RHESSI maps at 50-100 keV and 100-200 keV (Fig. |2J top right and bottom left) 
show double sources located at regions of opposite magnetic polarity (Fig. [2J top right) with 
one of the sources lying at the edge of the bright region in 195 A (Fig. El top left). Therefore, 
the HXRs in these energy ranges must be emitted near the footpoints of the loop. Note that 
the eastern HXR source is brighter than the western one. The magnetic field strength at 
the eastern source (~ +410 Gauss) is weaker than that of the western one (~ —510 Gauss). 
Since a stronger HXR source indicates a more efficient electron precipitation, the spatial 
relationship between the HXR sources and the magnetic field str ength can be interpreted as 



the result of magnetic mirroring of the HXR emitting electrons ( )Sakaolll994l ). 



The microwave source is located between the footpoint HXR sources (Fig. [21 bottom 
left). It is close to the coronal HXR (10-20 keV) source. The microwave source is also located 
at the region with weaker degree of polarization (Fig. [21 bottom right). We confirm from 
the spectral analysis (§ 12 . 3j) that the microwave emissions above 17 GHz are optically-thin 
non-thermal gyrosynchrotron emissions. Therefore, both the polarization information and 
the configuration of the longitudinal magnetic field indicate that the magnetic field at the 
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microwave source is quasi-perpendicular to the line of sight, which thus corresponds to the 
loop-top. Such a spatial distribution of the HXR and microwave emissions can be explained 
by the TPP model, if the microwaves are emitted by electrons trapped in the loop-top and 
the HXRs are emitted by electrons precipitating into the footpoints. 



2.3. Spectra 

We analyzed the temporally-resolved (but spatially-unresolved) spectra during the im- 
pulsive phase of the flare. We fit the RHESSI 40-250 keV spectrum at each time interval 
with a double power-law function of the form 

fUK-l <t)^ {t) , ife<eb(f), 

b(t)e-T*® ife>e b (t), [L) 

where e is the photon energy, 7 L (t) and 7 H (t) are the spectral indices of the lower- and 
higher-energy parts, e\,{t) is the break energy, and b(t) = a(t)eb(t) 7H ^ _TL ®, respectively. 
The upper panel of Figure [3] shows an example RHESSI energy spectrum and its fitting 
result. We chose energy bins of 2 keV from 40 to 60 keV, 2.5 keV from 60 to 100 keV, 
5 keV from 100 to 150 keV, and 12.5 keV beyond 150 keV, and a temporal resolution of 
4 s (approximately equal to the RHESSI rotation period). We used the front segments of 



detec tors #3, 4, and 8, which have the best energy resolution below ~ 100 keV (ISmith et al. 



20021 ). For convenience of analysis, the range of (t) was limited to 70 to 130 keV. 



We also obtained the temporally-resolved microwave spectral index from the NoRP 
data. After integrating it by 2 s to improve the statistics, we fit the NoRP spectrum taken 
with five frequencies of 2, 3.75, 9.4, 17, and 35 GHz at each time interval with a generic 



function (ISilva et al.ll2000l ) 



g{u,t) = a x {t)v a ^ [l-exp(-a 3 (tK a4(<) )] (2) 

ai(tK 2W , if ^< turnover (*), 

ai (t)a 3 (t)lS-^- a ^\ if V > IWoverW, 

where v is the frequency. We obtained the best-fit spectral index of the microwave flux 
density in the optically thin (higher frequency) regime, a(t) = a±(t) — a 2 (t) (positive value), 
as well as the turnover frequency ^turnover (0- The lower panel of Figure [3] shows an example 
NoRP microwave spectrum and its fitting result. We confirm that turnover (0 is l ess than 17 
GHz during the impulsive phase, indicating that the microwave emissions above 17 GHz are 
certainly optically thin, non-thermal gyrosynchrotron emission. 

Figure H] (upper) shows the time profiles of the spectral indices of the non-thermal 
emissions. The blue and red asterisks are the spectral indices of the lower-energy ( < 100 
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keV) and higher-energy ( > 100 keV) HXRs (hereafter 7£ bs (t) and 7H bs (£)), an d the green 
diamonds are the spectral indices of the microwaves in the optically thin regime (hereafter 
a ohs (t)), respectively 

The time profile of 7£ bs (i) shows the so-called soft-hard-soft behavior for each spike 
except the last. For example, 7£ bs (t) is ~ 4.8 at 01:02:10 UT, becomes hard (~ 3.7) at 
01:02:40 UT, and softens again (~ 4.0) at 01:03:10 UT. In addition, 7£ bs (t) also shows the 
so-called soft-hard-harder behavior during the entire course of the impulsive phase. However, 
the time profile of 7H bs (t) behaves differently from that of 7£ bs (t). Neither the soft-hard-soft 
behavior nor the soft-hard-harder behavior can be seen in the time profile of 7H bs (t). The 
values of the microwave spectral index are quite smaller than those of the HXR spectral 
indices. This implies that the inferred energy distribution of the microwave emitting electrons 
is harder than that of the HXR emitting ones. We will carry out a spectral analysis in § H] 
to confirm whether the energy distribution of the microwave emitting electrons is actually 
harder than that of the HXR emitting ones. 

In addition to this, we find that the time profile of 7H bs (t) is similar to that of a ohs (t), 
although the absolute values of their spectral indices differ by ~ 2. We also find that the 
time profile of 7H bs (^) (and a obs (t)) is delayed from that of 7£ bs (£). This tendency is especially 
seen during 01:03 - 01:05 UT. The cross-correlation functions of the spline-interpolated time 
profiles of the spectral indices are shown in the lower panels of Figure HI We find: (1) the 
time profile of a ohs (t) and that of 7H bs (^) show a peak correlation without a time delay (within 
a temporal resolution of 4 s); and (2) the time profile of 7H bs (^) is delayed by about 10 s from 
that of 7£ bs (£). The similarity of the time profiles between 7H bs (t) and a ohs (t) indicates that 
the higher-energy HXR emitting electrons and the microwave emitting electrons are from 
the same population and in a similar energy range. The delay of the time profile of 7H bs (^) 
from that of 7° bs (t) ma y be interpreted as an electron energy- dependent time delay (e.g., 



Aschwanden et al.lll997l ). 



3. A model for electron transport 



We modeled the observed spectral behavior in the 2003 May 29 flare in terms of TPP. 



The TPP model has been implemented in various 


ways, analytic (Me 


rose & Brownl 1976: 


Mad 


Gnnon et al.lll983: MacKinnon 19861. 


19881. 1991 


:lLeach & Petrosian 


1981 


: Lu & Petrosian 


1988 


and numerical approaches ( 


McClementd 


1990al 


bl:lHamilton et al. 199ol: I 


dacKinnon & Craie 


1991 


Fletcher & Martens 1998: 


Lee & Gary 


2000) to treat electron tr 


ansport in the flare 
vs an explicit treat- 


loop. We follow the approach made by 


McClementsI dl990b) that al\o\ 



ment of electron pitch-angle diffusion in a spatially-homogeneous magnetic loop. 
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3.1. Basic equations 



We trea t the sp atially-homogeneous, gyro-averaged Fokker-Planck equation following 
McClementsI Jl9901 



ON d 



+ 



dt dE 



EN 



F(E,fi,t) = j- 



D 



/</< 



dN 
<9/i 



Q{E,v,t). 



(3) 



Here, N(E,/j,,t) is the trapped electron distribution (number of electrons per unit energy 
per unit pitch-angle cosine), Q(E,fi,t) is the electron flux (number of electrons per unit 
energy per unit pitch-angle cosine per unit time) injected into the trap region, F(E,fi,t) is 
the electron flux precipitating into the footpoints, E and D m are the Coulomb energy loss 
rate and pitch- angle diffusion coefficient, E = T — 1 is the kinetic energy in units of the 
electron rest mass energy m e c 2 , T is the Lorentz factor, m e is the electron mass, c is the 
speed of light, and /i is the pitch-angle cosine, respectively. 

We adopt the Coul omb energy loss rat e and pitch-angle diffusion coefficient for a fully 
ionized plasma given by iLeach fc Petrosianl (1198 IT ) , 



E 
D 



-Kn/P 
Kn , 



(3 3 T 2 



-u E E, 



H 



(K = Ancrl In A), 

-/0> 



(E + 2) 



(4) 
(5) 



where r = 2.82 x 10 -13 cm is the classical electron radius, n is the ambient plasma number 
density, In A ~ 25 is the Coulomb logarithm for the typical solar coronal condition, z/g = 
Kn/((3E) is the Coulomb collision frequency, and (3 = a/1 — T~ 2 . The ambient plasma 
number density and the Coulomb logarithm are treated as constant in this paper. We 
neglect other Coulomb diffusio n coefficients such as Dee and D^e which are smaller than 



Dfj,u by a factor of order In A (IHamilton et al.lll990l ). Here, any other physics of electron 



kinematics such as wave-particle interactions are ignored for simplicity. 

For the precipitating electron flux F(E,fi,t), we again follow iMcClements (|l990aB ). 



F{E,n,t) 



L/Mv{E)\ 
H{\n\ - He) 



N{E,v,t) 
N{E,fjL,t), 



(6) 



where L is the characteristic scale length (i.e., loop length), H is the Heaviside step function, 

1 /2 

He is the loss cone angle cosine, v(E) = c(3 = c [l — (1 + E)~ ] is the velocity of an 
individual electron, and t c (E,h) is the electron loop-transit time along a magnetic field 
line, respectively. Here we assume a symmetrical magnetic loop with abrupt increases of 
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magnetic field intensity and ambient plasma number density below the mirror points near 
the footpoints, and we thus assume symmetrical precipitation in each h alf of ix space (n < 
and yU > 0). This form of equation is for mathematical convenience ( [MacKinnon! Il988l ). 



We take the parameters as follows: L 



3 x 10 9 



cm, n = 3 x 10 



10 



cm 



and /i c = 0.7. 



We adopt the value for the characteristic scale length L based on the distance of the HXR 
sources in the 2003 May 29 flare, assuming a semi-circular shape for the loop. The number 
density n in our model is slightly smaller than the observed value of (6 — 8) x 10 10 cm -3 
derived from the GOES observation during 01:03 - 01:05 UT (assuming a volume of L 3 and 
a filling factor of unity). We assume that electrons are trapped in an outer loop with a 
lower density than the brightest SXR loop. In the outer loop, energy dissipation (electron 
bombardment) at the footpoints and the resultant filling with evaporated chromospheric 
plasma have not yet occurred, whereas they have already occurred in the inner SXR loop. 
The value for the loss cone angle cosine /i c corresponds to a magnetic mirror ratio of 2, 
that is, the ratio of magnetic field strength at the footpoint to that at the trap region. We 
ad opt this value based on the statis tical analysis of flare data taken with CGRO and Yohkoh 
by lAschwanden et al.l (119981 . Il999l ). They estimated the fraction of directly-precipitating 
electrons to trap-precipitating ones and derived a magnetic mirror ratio of 1.2 — 3. 



3.2. Time-dependent injection flux 



We give the time-dependent, single power-law electron injection flux Q(E,/i,t), 



Q{E,fM,t)=A(t)[^- 



-<5i„(t) 



0(/x), 



(7) 



where we adopt a pivot point energy (IGrigis fc Benzl 120061 1 Eq of 0.098 (= 50 keV), and 



A(t) 
5 in (t) 



exp 



t-30 
25 



exp 



t-90 
25 



, Tit 

4.5 + cos^ ( — 
60 



(9) 



for < t < 120 sec. Equations (|HJ) and (Q describe double peaks electron injection with 
a soft-hard-soft spectrum. We adopt this form for the electron injection because the HXR 
spectrum in the lower-energy regime observed in the 2003 May 29 flare shows the soft-hard- 
soft behavior. 

The remaining term of 0(/x) gives the (time- and energy-independent) pitch-angle dis- 
tribution of the injection flux. We perform calculations for two cases of the pitch-angle 
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distribution: the pancake distribution, <fi(ji) = exp[— (/^/0.5) 2 ], and the isotropic distribu- 
tion, = const. 

Using equations PJ - (JSJ, we numerically solve equation ([3]). We apply a finite difference 
method with operator splitting. The differential operators in equation ([3]) are split into two 
terms: the diffusion term and the remaining terms. We solve the diffusion term by using the 
Crank-Nicholson method with central difference, accurate to second order in time. We em- 
ploy the symmetric boundary condition in [i space which satisfies total number conservation. 
We set 97 grid points in fi space by the following manner, 

Ho = 0, Afi = 0.016, 

fi n = + A/i(l - in = 1,2, • • • ,95), (10) 

^96 = 1, 

giving coarse grids for smaller \x and fine grids around the loss cone. For th e remaining term 



we use an analytic solution given by the method of characteristics (see, e.g.. lCraig et al.lll985 



MacKinnonlll986l ). For the necessary interpolation at the intermediate location between the 
grid points in energy space, we use a single power-law function. We set 256 grid points in 
energy space, logarithmically-spaced from 50 to 5000 keV. 



4. Calculation results and discussion 

In this section, we present our calculation results of the Fokker-Planck equation. Figures 
[5] and [6] show the electron distribution in phase space (E,/a), calculated for the pancake and 
isotropic pitch-angle distributions of the injection, respectively. The left panels show the 
trapped electron distribution N(E,fi,t) at selected times of t = 10 sec (top) and t = 50 sec 
(bottom). The right panels show the slope s of N(E,^,t) in energy at the selected times, 
derived from the following equation, 

log [N(E l+1 ,H,t)/N(E,fi,t)} 

log(iW^) ' 1 J 

where subscript % denotes the grid position in energy space. We can see from these right 
panels that the electron energy distribution outside the loss cone (\x < 0.7) is harder in the 
decay phase (t = 50 sec) than in the rise phase (t = 10 sec). The electron energy distribution 
inside the loss cone (/i > 0.7), which is related to the precipitating electron flux F(E,fi,t) 
by equation ([6]), is softer than that outside the loss cone. The electron distribution inside 
the loss cone shows different features between the pancake and isotropic cases. 



In § 14.11 we discuss the time evolution of the trapped electron distribution and the 
precipitating electron flux in energy space by using the double power-law fitting method. 
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Our calculations confirm the electron trap and precipitation regardless of the weak and 
strong diffusion limits. In § 14.21 we discuss the time evolution of the non-thermal emissions 
predicted by the electron model for comparison with the observations. 



4.1. Evolution of electrons 

To illustrate the evolution of the trapped and precipitating electrons, we introduce the 
pitch-angle integrated variables: 



N„(E,t) = [ diiN(E,fi,t), 
Jo 



F„(E,t) = I dfiF(E,fi,t)= I rf/x ^f'^ . (12) 



Next, we fit the electron energy spectra from 50 to 400 keV at each time by the double 
power-law function of equation (CD) (but replacing [e, 7] by [E, A] for N^E.t), and by [E, 5} 
for F^Ejt)). The range of E^(t) is limited to 100 to 200 keV for convenience. An example 
of our fitting of the electron energy spectra is shown in Figure [71 In this way, we obtain the 
spectral indices Al(£), An(t), ^l(^), and 5u(t). 

Figure M shows the time profiles of the spectral indices of the electrons in the pancake 
case. The left panel shows Al(£) (lower-energy regime, solid line) and An(t) (higher-energy 
regime, dashed line). The break energy is ~ 170 keV. The right panel shows S^(t) (lower- 
energy regime, solid line) and 5n(t) (higher-energy regime, dashed line). The break energy 
is ~ 130 keV. 

The spectra of F^E^t) are quite softer than those of N^{E,t). To understand this in 
terms of a diffusion regime, we introduce the precipitation rate 

^ = ^ oc E\ (13) 

In the strong and weak diffusion limits, the precipitation rates are respectively evaluated as 
(using the non-relativistic expression), 

l/r e oc v oc E - 5 , (strong) 

l/r d ~ oc E" 1 - 5 , (weak) ( ' 

that yields — 1.5 < x < 0.5. Whether diffusion for an electron is weak or strong depends on 
its energy E. This means that x itself is a function of E. When > r c (which corresponds 
to higher E), the diffusion is weak and x takes a negative value, and vice versa. In our 
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situation, the energy at which ra ~ r c is about 20 keV (using \i = 0.7). All of the electrons 
in our calculation (> 50 keV) satisfy t$ > r c , and so are in the weak diffusion regime. 
Because t& becomes much longer than r c for the higher-energy electrons, x approaches its 
lower limit (weak diffusion limit) of —1.5 with increasing E. We find that —1.5 < x < —1.0 
in the high-energy regime and x ~ —0.7 in the lower-energy regime from Figure This is 
consistent with the previous statement, and thus confirms our general treat ment of the TPP 



model in the weak diffusion regime considered by lMelrose fc Brown! ( 119761 ). 



In the left panel of FigureEl the time profile of Al(£) shows the soft-hard-soft behavior in 
each spike (0 < t < 60 and 60 < t < 120), achieving its hardest values around t = 30 and t = 
90. However, the time profile of A H (i) does not show such behavior but shows the soft-hard- 
harder behavior in each spike. The energy-dependent trap efficiency yields this difference 
of the temporal variation of the spectral indices between the lower and higher energies. 
The trapped electrons are lost via Coulomb energy loss and precipitation. The Coulomb 
energy loss rate is smaller for the higher-energy electrons. As mentioned above, the weak 
diffusion yields a precipitation rate that is also smaller for the higher-energy electrons. This 
means that the escape time scale from the phase space becomes longer for the higher-energy 
electrons. For the lower-energy electrons, the escape time scale is on the order of 1 sec. This 
is much shorter than the injection time scale (see eq. ([8])). Therefore, Al(£) reflects the 
temporal variation of the spectral index of the injection flux 5i n (t) described by equation 
(Q. The escape time scale becomes comparable to or longer than the injection time scale for 
the higher-energy electrons. The higher-energy electrons stay at the trap region and their 
energy remains high for a longer time. As a result, the spectrum in the higher-energy regime 
becomes harder in the decay phase of each spike (30 < t < 60 and 90 < t < 120). 

Figure [9] shows the time profiles of the spectral indices of the electrons in the isotropic 
case. The left panel shows Al(£) (lower-energy regime, solid line) and An(t) (higher-energy 
regime, dashed line). The break energy is ~ 170 keV. The right panel shows S^(t) (lower- 
energy regime, solid line) and 5u(t) (higher-energy regime, dashed line). The break energy 
is ~ 120 keV. 

There is no much difference between the left panels of Figures [H] and O indicating that 
the evolution of the trapped electron distribution in energy space is almost independent 
of the pancake or isotropic pitch-angle distribution of the injection. We find a difference 
between the right panels of Figures [8] and [91 In both cases, the time profiles of S^(t) show 
the soft-hard-soft behavior in each spike. In the isotropic case (right panel of Fig. [9]), the 
time profile of Sn(t) shows the soft-hard- soft behavior same as that of S^it). In the pancake 
case (right panel of Fig. [8]), however, the time profile of Suit) is delayed from that of S^it). 
We interpret these features in terms of the difference of the injection pitch-angle distribution. 
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The precipitating elect rons consist of two diffe rent types: directly-precipitating and 



trap-precipitating electrons (lAschwanden et al.lll998l ). Electrons injected with small pitch- 



angle directly precipitate without being trapped, while those injected with large pitch-angle 
are trapped once and subsequently precipitate via pitch-angle scattering. When the injected 
electrons have an isotropic pitch-angle distribution, the precipitating electrons include the 
directly-precipitating ones as well as the trap-precipitating ones. When the injected electrons 
have a pancake pitch-angle distribution, on the other hand, almost all of the precipitating 
electrons are trap-precipitating ones. The right panels of Figure which show the slope 
of N(E, fj,, t) in energy in the pancake case, show a "propagating feature" (contours around 
the loss cone) from outside the loss cone to inside the loss cone. This corresponds to the 
trap-precipitating electrons. The "propagation speed" is faster for the lower-energy electrons 
because it is governed by the pitch-angle diffusion coefficient (eq. (J3J). This feature is less 
clear in the isotropic case (right panels of Fig. This is because in the isotropic case the 
trap-precipitating electrons merge with the directly-precipitating ones. 

The directly-precipitating electrons precipitate on a time scale of r c . The precipitation 
time scale of the trap-precipitating electrons, ~ ra, is longer than r e . A fraction of these two 
components, which is determined by the pitch-angle distribution of the injection, determines 
the precipitation time scale of all of the precipitating electrons. The precipitation time scale is 
longer in the pancake case than in the isotropic case. For the lower-energy electrons, however, 
the precipitation time scale is shorter than the injection time scale in both the pancake and 
isotropic cases. Therefore, S^it) reflects the temporal variation of S in (t) regardless of the 
injection pitch-angle distribution. For the higher-energy electrons, the precipitation time 
scale becomes longer due to the weak diffusion. It is comparable to or longer than the 
injection time scale in the pancake case, but not in the isotropic case. Consequently, the 
time profile of 5h(^) is delayed from that of Si n (t) (and thus S^(t)) only in the pancake case. 



4.2. Evolution of radiation 



In this section, we show the time evolution of the non-thermal emissions predicted by 
the Fokker-Planck calculation results, and compare them with the observations. The spa- 
tial distribution of the non-thermal emissions in the 2003 May 29 flare shows a loop-top 
microwave source and double footpoint HXR sources (Fig. [2]). This supports the inter- 
pretat ion that the trapped elect rons N(E,fi,t) emit microwaves via gyrosynchrotron radi- 
ation (|Ramatyill969l ; lDulk!ll985f) and the precipitating electrons F(E,fi,t) emit HXRs via 
thick-target bremsstrahlung (|Brownlll97ll ). We numerically calculate the thick-target HXR 
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intensity at a photon energy e by 

1 



I(e,t) 



AtxR 2 



dEF^(E, t) 



E dE ,nv(E')a B (e, E>) 



E'ue(E') 



(15) 



where 0"B(e, E) is the direction- integrated bremsstrahlung cross section given by lHaugl (119971 ). 
ue is given in equation (jlj), and R = 1 AU. We use the pitch-angle integrated electron flux 
because electrons precipitating into the thick-target region are quickly isotropized. 

For the microwave emission, we numerically calculate only the gyrosynchrotron emissiv- 
ity from the t r apped electrons N(E, fi, t) by using an approximate analytic expression given 
by iPetrosianl (119811 ) (see Appendix). Although a general description of gyrosynchrotron 
radiation includes absorption, this approximation is sufficient for our purpose because we 
discuss the microwave spectral behavior only in the optically-thin regime. We calculate the 
gyrosynchrotron emissivity in a harmonic range of 10 — 100 with a nominal viewing angle 
= 75°. When the magnetic field strength at the emission site is on the order of 100 Gauss, 
the frequency range in our calculation corresponds to the ranges of NoRP and NoRH. 



4-2.1. Lightcurves 

Figure [10] shows the lightcurves of the non-thermal emissions in the pancake case. The 
upper part shows the HXR fluxes at 50 keV, 99 keV, and 197 keV on a normalized scale. 
The dashed lines denote the peak times of each spike in the HXR 50 keV flux. We can see 
that the peaks of the HXR lightcurves within 50-200 keV are almost coincident. The lower 
part of Figure [10] shows the microwave emissivities at 17 GHz and 35 GHz on a normalized 
scale, assuming a magnetic field strength B at the trap region of 300 Gauss. The peaks of 
the microwave emissivities are delayed from the HXR flux at 50 keV by about 5 s. These 
tendencies are also observed in the isotropic case, and are consistent with the non-thermal 
lightcurves for the 2003 May 29 flare (Fig. [JJ). 



4-2.2. Spectra 

To illustrate the spectral variation of the calculated emissions, we fit the calculated HXR 
spectrum within 50-200 keV at each time with the double power-law function of equation 
([TJ). The range of 6b(t) is limited to 75 to 125 keV. We also fit the calculated microwave 
spectrum from 17 to 35 GHz at each time with a single power-law function. 

Figures [UJ and [12] show the time profiles of the spectral indices of the non-thermal 
emissions. In these figures, the upper plots are the calculation results in the pancake and 
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isotropic cases, respectively. The lower plots are the observation results (but the microwave 
spectral indices are multiplied by a factor of 2) during 01:03 - 01:05 UT for comparison 
with the calculations. The complete set of the observation results is in Figure HI Colors 
(blue, red, and green) denote the lower-energy HXRs, the higher-energy HXRs, and the 
microwaves, respectively. Hereafter, spectral indices of the lower-energy HXRs, the higher- 
energy HXRs, and the microwaves in the calculated spectrum are named 7 FP (^)> 7h P (^); an d 
a FF (t), respectively. 

Our model calculation results show some agreements with the observations. The values 
of the microwave spectral indices are smaller by ~ 1.5 — 2 than those of the HXR spectral 
indices in both the pancake and isotropic cases. This result is quantitatively consisten t 
with not only our observations but also with previous reports (e.g., ISilva et al.lll997l . |2000| ). 



This is understood by using simple analytic formulae that relate the spectral index of the 
emissions to that of the parent electrons. We assume tha t the trapped electrons have a 
power-law energy distribution, N^E) oc E~ A . iDulkl fll985h gives an empirical relationship 
of the spectral indices between the gyrosynchrotron emission in the optically thin regime 
and the parent electron energy distribution at the site, 



a 



0.9A - 1.22. 



(16) 



Although this relationship is derived under the assumption that the parent electrons have an 

isotropic pitch-angle distribution, we use it for a rough estimate of the microwave spectral 

index from electrons with an arbitrary pitch-angle distribution. For HXRs, we can use 

equation ([TBI that shows the energy distribution of the precipitating electron flux, F^(E) oc 
£-(A-*)_ An analytic 

expression for non-relativistic thick-target bremsstrahlung gives the 
relationship of the spectral indices between th e HXR emissio n and the parent electron flux 
precipitating into the thick-target region (e.g., 



Hudson! 11972( 1. 



7 = (A -x) - 1. 



;i7) 



Subtracting equation ([16]) from equation (TPTj) . we find the difference of the spectral indices 
between the HXR and microwave, 



7 



a 



O.lA-x + 0.22. 



In the weak diffusion regime x takes negative value, around —1.0 in our calculation. This 
eventually yields 7 — a ~ 1.5 (for A = 3) which is in agreement with the observations. Thus 
we conclude that the difference of the spectral indices between the HXR and microwave 
emissions can be interpreted as a consequence of the parent electron transport in the TPP 
model in weak diffusion regime. 
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There is no much difference in the time profiles of a FP (t) between the pancake and 
isotropic cases, in that both show the soft-hard-harder behavior during the first spike (0- 
60 sec). Our calculations of the TPP model can successful ly reproduce the often observed 



soft-hard-harder behavior of the microwave spectrum (e.g.. ISilva et al.ll2000l ). On the other 
hand, there is a difference in the time profiles of 7l P (£) and 7h P (£) between the pancake and 
isotropic cases. In the isotropic case, the time profiles of 7n P (t) and 7l P (^) both show the 
soft-hard-soft behavior. In the pancake case, however, 7h P (£) shows similarity with a FP (t) 
rather than with 7l P (£), which only shows the soft-hard-soft behavior. During the decay 
phase of the earlier spike in the observation (01:03:30 - 01:04:00 UT), both 7 H bs (t) and a ohs (t) 
show hardening whereas 7£ bs (i) shows slight softening. During the rise phase of the later 
spike in the observation (01:04:00 - 01:04:20 UT), ^ s (t) and a ohs (t) show softening. These 
tendencies can be seen during the respective 30-60 sec and 60-80 sec only in the pancake case 
calculation (Fig. [TT]) . These qualitative agreements suggest that it is reasonable to consider 
the pancake pitch-angle distribution of the injection flux rather than the isotropic one to 
explain the observed spectral behavior in the 2003 May 29 flare. 



4.3. Validity of our interpretation 

In our current study, we have utilized the model electron distributions calculated for 
only two different types of the injection pitch-angle distribution (pancake and isotropic) with 
nominal values of parameters: the ambient plasma number density n and the loss cone angle 
cosine fi c (the mirror ratio). Since these parameters also affect the evolution of electrons, 
we have to systematically investigate the set of parameters with which we can reproduce 
the observation of the 2003 May 29 flare. We believe that n of an o rder of 10 10 cm~ c 



from the GOES observation and the mirror ratio of 1.2 — 3 derived by lAschwanden et al. 



( 119981 ) are reasonable ranges. A much higher mirror ratio drastically reduces the number of 
precipitating electrons and the resultant thick-target HXR emissions at the footpoints, which 
may be in disagreement with the observations. A much higher density (n > 10 11 cm -3 ) at 
the trap region would produce strong coronal thin-target HXR emissions, which is rarely 
observed. As such, we consider that a relatively lower density at the trap region (less than 
that of the SXR bright loop) and a relatively smaller mirror ratio should be reasonable, and 
thus we believe the values adopted in our calculation, n = 3 x 10 10 cm" 3 and /x c = 0.7, are 
in reasonable ranges. Our interpretation made earlier in this paper remains valid for the 
specified ranges of the number density and the mirror ratio. To achieve a better agreement 
between the calculation and the observation, we need to further refine the model distribution 
of the injection flux. 
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Lee et al.l (120001 ) and lLee fc Gary! (120001 ) reported the microwave observation of a GOES 
C2.8 flare on 1993 June 3. They further performed the calculations of the TPP model. They 
carried out a systematic investigation by varying the number density and injection time scale 
as well as the injection pitch-angle distribution, to search for the best parameter set that 
agreed with their observation. They concluded that the electrons are confined to a narrow 
range < 0.26) of pitch-angle and are injected into a low density (n~4x 10 9 cm -3 ) trap 
region. This number density is much lower than that we assumed in this paper. This may 
be because we observed the large X-class flare while they observed the small C-class flare. 

Our interpretation of the difference of the spectral indices between the HXRs and mi- 
crowaves based on the TPP model is subject to the observation that microwaves are emitted 
at the loop-top whereas HXRs are emitted at the footpoints. Our assumption in the TPP 
model that the HXR and microwave emitting electrons are treated separately is inadequate 
for a flare which shows, for example, microwaves as well as HXRs at the footpoints. In such a 
flare, the microwave emitting electrons are identical to the HXR emitting ones. This means 
x = 0.5, yielding a ~ 7 in equation (fl8|) . The footpoint microwave emissio ns would thus be 



expec ted to have almost the same spectral index as the footpoint HXRs. lYokoyama et al. 



( 120021 ) reported that the footpoint microwave emission has a softer spectral index than the 
loop-top one by ~ 2. Their result is consistent with our interpretation. Simultaneous obser- 
vations of HXRs and microwaves in a flare showing footpoint HXR and microwave emissions 
would be useful to further explore the validity of our interpretation. 



5. Conclusion 

We presented the comparative study of the non-thermal emissions of the flare occurring 
on 2003 May 29 using the RHESSI HXR and Nobeyama microwave observations. Further, 
we considered the electron transport model, TPP, to explain the observations. 

The 2003 May 29 flare showed two non-thermal HXR sources at the footpoints and a 
microwave source at the loop-top, as observed with RHESSI and NoRH. We interpreted this 
in terms of the TPP model. We presented the time profiles of the spectral indices of the 
higher-energy HXRs 7H bs (^) as well as the lower-energy HXRs 7£ bs (t) and microwaves a ohs (t). 
The spectra of microwaves and HXRs imply that the microwave emitting electrons have a 
harder energy distribution than the HXR emitting ones. We found that the time profile of 
7H bs (£) shows similarity with that of a ohs (t) rather than with 7£ bs (£), and is delayed from 
that of 7£ bs (t). 

We numerically solved the spatially-homogeneous Fokker-Planck equation for the TPP 



-18- 



model to describe the evolution of electrons. Precipitating electrons have a softer energy 
distribution than the trapped ones in the weak diffusion regime. Differences of the injection 
pitch-angle distribution especially affect the evolution of the precipitating electrons. 

We calculated the microwave and HXR emissions from the calculated trapped electron 
distribution and precipitation flux for comparison with the observations. The TPP model in 
the weak diffusion regime can yield a soft HXR spectrum and a hard microwave spectrum. 
The calculated difference of the spectral indices between the HXRs and microwaves, ~ 
1.5, is in agreement with the observations. We further found that a pancake pitch-angle 
distribution for the injected electrons rather than an isotropic distribution is more adequate 
to qualitatively explain the temporal variation of 7n bs (^)- By comparing the model calculation 
with the observation, we can constrain the pitch-angle distribution of the injected electrons, 
which is crucially important for understanding the electron acceleration mechanism in solar 
flares. 

Currently, we are improving our treatment of the TPP model to include the spatial 
inhomogeneity in the Fokker-Planck equation. Using this, a systematic investigation of the 
best parameter set to explain the observation is in progress, and will be reported in the 
future. 
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A. Calculation of the gyrosynchrotron emissivity 



A general calculation of the gyrosynchrotron radiation in a magnetized plasma (IRamaty 

19691 ) includes effects such as self-absorption, absorption by ambient plasma, and Razin 



suppression. These effects significantly contribute at low harmonics (v/ub ^ 10, where ub is 
the electron gyrofrequency) of the gyrosynchrotron radiation. In our current study, however, 
only the optically-thin gyrosynchrotron radiation at high harmonics (10 < vjv-Q ^ 100) from 
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mildly relativistic electrons (r < 10) is of i nteres t. Under such limited conditions, there 
is an useful expression given by iPetrosianl (|198ll ). We adopt his formula to predict the 
microwave emission from the Fokker-Planck calculation results. The approximate expression 
of the gyrosynchrotron emissivity at a frequency v and a viewing angle 9 with respect to 
the magnetic field, from mildly relativistic electrons with arbitrary energy and pitch-angle 
distributions, is as follows; 



e 2 v B 



v 



% ub sin 2 9 

where e is the elementary charge, and 



)POO pi 
J dT J d/iN(T -l,/i,t)Y(9,T,/j,)Z 2m (9,T,ij<), (Al) 



Y 



Z 



{cos 9 - ft/3) 2 + (1 - z 2 )(l - /3/xcosi 

(1- £2)1/2(1-/^ C os fl) 

zexp [(1 -z 2 ) 1 ' 2 ] PsmG(l-i 
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Fig. 1. — Lightcurves for the 2003 May 29 flare. Upper. The GOES SXR lightcurves during 
00:50 - 01:50 UT. The period during 01:00 - 01:05 UT (dot-dashed lines) is defined as the 
impulsive phase. Lower. Lightcurves of the non-thermal emissions during the impulsive 
phase. The upper part shows the HXR lightcurves (in units of counts/cm 2 /sec/keV) taken 
with RHESSI in three energy bands: 50-70 keV, 70-100 keV, and 100-200 keV from top to 
bottom. The lower part shows the microwave lightcurves (S.F.U.) observed in the NoRP 17 
GHz and 35 GHz bands from top to bottom. The dashed lines denote the peak times of each 
spike in the 50-70 keV band. 
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Fig. 2. — Spatial distribution of the emissions during the impulsive phase of the 2003 May 
29 flare (around 01:04:30 UT). The top left map is the TRACE 195 A image overlaid by the 
RHESSI 10-20 keV contours. The top right map is the SOHO/MDI magnetogram (in units 
of Gauss) overlaid by the RHESSI 50-100 keV contours. The dashed lines denote magnetic 
neutral lines. The bottom left map is the NoRH 34 GHz brightness temperature (in units 
of kelvin) image overlaid by the RHESSI 100-200 keV contours. The RHESSI images are 
reconstructed using the PIXON algorithm. Contour levels are 10%, 20%, 30%, 50%, 70%, 
and 90% of the peak intensity in each image. Thick and dashed contours overlaid on the 
NoRH 34 GHz brightness temperature image in the bottom right map are degree of right- 
and left-circular polarization at 17 GHz, respectively. Contour levels are 10%, 30%, 50%, 
and 70%. The angular resolution (beam size) of the NoRH 34 GHz map is ~ 9". The dotted 
lines denote heliographic grids in 2° increments. 
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Fig. 3. — HXR and microwave spectra taken with RHESSI and NoRP during the 2003 May 
29 flare. The upper plot shows the HXR photon spectrum (points with error bars) and the 
fitted double power-law function (solid line) during 01:04:28 - 01:04:32 UT. Values of 7l, 7h> 
and e b (keV) are determined to be {3.01 ± 0.001, 3.88 ± 0.003, 93.0 ± 0.1}. The lower plot 
shows the microwave spectrum at five frequencies of 2, 3.75, 9.4, 17, and 35 GHz (asterisks) 
and the fitted model described by equation (J3J) (solid line) at 01:04:29 UT. Values of a and 
J-Wiovcr (GHz) are determined to be {1.25, 11.1} 
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Fig. 4. — Upper. Time profiles of the spectral indices of the non-thermal emissions during 
the impulsive phase of the 2003 May 29 flare. The blue, red, and green symbols are 7£ bs (t), 
7^ bs (t), and a ohs (t), respectively. The break energy is ~ 100 keV. The solid line denotes the 
RHESSI count flux at 50-100 keV. Note that some data which have large uncertainty are 
omitted in the plot. Lower. Cross-correlation functions of two time profiles of the spectral 
indices as a function of time lag, during 01:03 - 01:04 UT (left) and 01:04 - 01:05 UT (right). 
The solid lines are cross-correlation of 7£ bs (t) and ^^(t), an d the dashed lines are that of 
7g bs (£) and a ohs (t), respectively. 
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Fig. 5. — Electron distribution in phase space (E,fj,), calculated for a pancake pitch-angle 
distribution of the injection flux. The top left panel shows the trapped electron distribution 
N(E,fi,t) at a time of t = 10 sec. Contour levels are 10" 6 , 10~ 5 , 1(T 4 , 10~ 3 , 1(T 2 , and KT 1 
of the maximum value. The top right image with contours shows the slope of N(E, /i, t) 
in energy, determined from the ratio of N(E,(j,,t) at two adjacent energy grids (eq. (ITT]) ). 
Contour levels are 3, 4, 5, and 6. The dotted lines denote the loss cone angle cosine. The 
bottom two panels are same as the top ones, but at a time of t = 50 sec. 
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Fig. 6. — Same as Figure [5l but calculated for an isotropic pitch-angle distribution of the 
injection flux. 
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Fig. 7. — Electron energy spectra in the pancake case. The upper panel shows the energy 
spectra of the trapped electron distribution N^E, t) at selected times of t = 10 sec (squares) 
and t = 50 sec (triangles) with the fitted double power-law function (solid lines). Values 
of A L , A H , and E h (keV) are determined to be {4.13 ± 0.017,4.51 ± 0.028, 165 ± 7.9} at 
* = 10 sec, and {3.73 ± 0.019,3.12 ± 0.025, 153 ± 4.7} at t = 50 sec. The lower panel 
shows the energy spectra of the precipitating electron flux F^{E,t). Values of 5^, 5u, and 
E h (keV) are determined to be {4.74 ± 0.023,5.52 ± 0.020, 133 ± 3.1} at t = 10 sec, and 
{4.45 ± 0.016, 4.46 ± 0.018, 100} at t = 50 sec. 
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Fig. 8. — Time profiles of the spectral indices of electrons in the pancake case. Left: Time 
profiles of Al(£) (lower-energy regime, solid line) and An(t) (higher-energy regime, dashed 
line). The break energy is ~ 170 keV. Right: Those of S^(t) (lower-energy regime, solid line) 
and 5n(t) (higher-energy regime, dashed line). The break energy is ~ 130 keV. 
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Fig. 9. — Same as Figure El but in the isotropic case. Break energies are ~ 170 keV (left) 
and ~ 120 keV (right). 
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Fig. 10. — Non-thermal emissions predicted by the Fokker-Planck calculation results in the 
pancake case. The upper part shows the HXR fluxes at 50 keV, 99 keV, and 197 keV on a 
normalized scale from top to bottom. The lower part shows the microwave emissivities at 17 
GHz and 35 GHz on a normalize scale from top to bottom. A magnetic field intensity of 300 
Gauss and a viewing angle 9 = 75° are assumed for the microwave emissivity calculation. 
The dashed lines denote the peak times of each spike in the HXR 50 keV flux. 
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Fig. 11. — Time profiles of the spectral indices of the non-thermal emissions. Upper. Calcu- 
lation result in the pancake case. The blue, red, and green lines with error bars are 7l P (^) ; 
7n P (t), and a FP (t), respectively. The break energy is ~ 90 keV. The solid line denotes the 
predicted HXR flux at 50 keV on a normalized scale. Lower. Observation result during 01:03 
- 01:05 UT. Note that the microwave spectral indices are multiplied by a factor of 2. For a 
complete set of these data, see Figure @J 
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